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Abstract 
> 

^" . We describe a general method that enables us to obtain all the singular terms of helicity 

amplitudes of n-parton processes at one loop. The algorithm uses helicity amplitudes 

at tree level and simple color algebra. We illustrate the method by calculating the 

^^ ■ singular part of the one loop helicity amplitudes of all 2 — > 3 parton subprocesses. The 

^. , results are used to derive the soft gluon limit of the cross sections of all 2 ^ 4 parton 

■^ I scattering subprocesses which provide a useful initial condition for the angular ordering 

o ■ approximation to coherent multiple soft gluon emission, incorporated in existing Monte 

JL . Carlo simulation programs. 
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1 Introduction 

The technical developments achieved recently in perturbative QCD calculations (helicity 
method fl], ||, H, Q, string theory based derivations [^, H) enable us to calculate one-loop 
corrections to helicity amplitudes up to five-parton and perhaps also to six-parton processes. 
Some partial results for the five parton amplitudes are by now available. 

In this development the utilization of the helicity method was decisively important. It was 
used both in the string theory based derivation of the one-loop corrections to the four-gluon 
[^] and five-gluon [Q helicity amplitudes, as well as in the recent conventional diagram- 
matic calculation of the one loop corrections of the helicity amplitudes for all 2 ^ 2 parton 
processes 0. 

The helicity method which was originally proposed for calculating tree amplitudes 0] 
can be extended to the calculation of loop corrections by observing that the difficulties given 
by the appearance of 75 in the definition of the helicity states can be circumvented if in- 
stead of conventional dimensional regularization the 't Hooft-Veltman regularization [|, ^ 



or dimensional reduction p!0| , pH] , |I2|, |T3| are utilized. While employing these regulariza- 
tions requires some care in calculating two- or more-loop contributions, their use is rather 
straightforward in the case of one-loop corrections. The regularization dependent terms in 
the loop contributions obtained with conventional dimensional regularization, dimensional 
reduction or 't Hooft-Veltman regularization are determined by the universal process inde- 
pendent structure of the singular contributions. The process independent transition rules 
have been worked out recently 0. 

An important technical issue in the evaluation of next-to-leading QCD corrections to 
physical cross-section is the analytic cancelation of the soft and collinear singularities between 
the loop and the Bremsstrahlung contributions. The singular terms in the cross-sections 
coming from Bremsstrahlung contributions have universal, process independent structure. 
They depend only on the color representation of the external legs and each leg contributes 
additively. It was recently shown that these terms can easily be obtained with the help of 



a Born-level algorithm |T^ (see also JT^). Since the virtual corrections are cancelled by the 
Bremsstrahlung contributions, one expects that a simple Born-level algorithm should also 
be sufficient to calculate directly the singular terms of the virtual corrections as well. The 
aim of this paper is to show that indeed such an algorithm existsQ. 

In section 2 we give a brief description of our algorithm for calculating the singular 
terms of one-loop helicity amplitudes. In section 3 we illustrate the power of the method by 
calculating the singular terms of one-loop helicity amplitudes of all 2 — »• 3 subprocesses. The 
soft singular terms of the one-loop helicity amplitudes of the 2 — >■ n — 2 parton processes 
determine also the soft limit of the cross sections of the 2 ^ n — 1 parton scattering processes 
via the cancelation theorem. In view of its phenomenological usefulness, we summarize our 
results for the soft limit of the Born cross sections of all 2 — i> 4 subprocesses in Section 4. 



^ A comprehensive treatment of the soft Umit in QCD is presented in|15|. For a recent paper in the large 



Nc hmit see |16 . 



2 Singular terms in one-loop helicity amplitudes 

In D-dimensional regularisation the virtual corrections to helicity amplitudes at one loop 
have soft 1/e^ and 1/e, coUinear 1/e and ultraviolet 1/e singularities. The ultraviolet singular 
term is trivially cancelled by the MS coupling constant counter term 

,UV in-2){ATiy f gY .tree 



where n denotes the number of the external legs, jSq is the first term in the /^-function 

/3o = (lliV,-2iV^)/3, (2.2) 

g denotes the renormalized coupling constant and D = A — 2e. Color and helicity labels are 
supressed. In the following we shall discuss only renormalized helicity amplitudes. 

In axial gauge the coUinear singularities come from the self energy corrections to the 
external lines [|T^. For each helicities and color subamplitudes they are proportional to the 
Born term since the Altarelli-Parisi functions, Pa/a{z), for diagonal splitting preserve helicity 
in the z — *> 1 limit. Therefore, the coUinear singularities have the form 

4T = -(£)'E^^"^^ (2.3) 

The constant 7(a) represents the contribution from virtual graphs to the Altarelli-Parisi 
kernel Pa/a{0 '^^ equivalently it is defined by the behavior of Pa/a{0 '^^ar ^ = 1 

£ d^ Pa/aiO = 2C(a) ln(l -z)+ 7(a) + O {1 - z) , (2.4) 

where C{a) is the color charge of parton a. Specifically, C{a) and 7(a) are 



C{g) = K , l{9) = iMj^, for gluons, 
*^(9) = ^kwr^ ^(^) = ^^^Sf^' fo^ quarks. 



(2.5) 



The structure of multiple soft emission from hard processes in QED was investigated by 
Gammer and Yennie |l^. They have shown that the energetic electrons participating in a 
hard process receive an eikonal phase factor. In quantum chromodynamics, the situation is 
very similar except that the eikonal factor is a matrix equal to the path order product of the 



matrix- valued gluon field [0, [T^, |2^ . 



For one soft gluon, the main result is very simple: it states that the singular contributions 
come from a configuration where the soft gluons are attached to external legs of the graphs. 
Therefore, the soft contribution can easily be calculated in terms of the Born amplitude. For 
example, let us consider the contribution of a virtual soft gluon to the matrix element of a 
process with n external legs, 2 — > (n — 2). 

The insertion of a soft gluon which connects the external legs i and j has a twofold effect: 
First, the Born amplitude gets rotated in the color space by the insertion of the color matrices 
appearing in the two vertices of the soft line connecting the hard lines i and j 

A{2^n- 2),,,,..e„ -^ Y. tlc't^cA^ ^n- 2),,..^,,, (2.6) 



where the labels of the initial partons are {n — 1) and n; q's denote color indices for the 
external partons and t" / is the SU(3) generator matrix for the color representation of line i. 
That is, t"j is (l/2)A°j- for an outgoing quark, — (l/2)A*j" for an outgoing antiquark, and ifaij 
for an outgoing gluon. For an incoming parton, one can use the same formula as long as one 
uses the conjugate color representations, —{1/2)X*°- for a quark, (l/2)A°j- for an antiquark, 
and ifaij for a gluon. 

Secondly, after carrying out the loop integral and dropping singular terms corresponding 
to coUinear configurations, we pick up the same eikonal factor as in QED, i.e. 

where we introduced the short hand notation 

Using unitarity this form of the soft factor can be confirmed by integrating the gluon mo- 



menta over the Bremsstrahlung eikonal factor []T4 . 

To obtain the complete soft singular contributions, one should consider all configura- 
tions where a gluon connects two external legs. So the soft singular terms of loop helicity 
amplitudes can be written as 

^ir(2 -. n - 2).,.,...„ = - E (£)' ^r^ (-7^) ^ ^ ^.S.^^-^'"'^ - - " 2)....e..e;..e„ 

(2.9) 
The color matrix of the Born term is an invariant tensor, therefore the sum of insertions to all 
external lines must vanish (soft identity [Q)- The soft structure given by ( |2.9| ) can easily be 



decomposed into color subamplitudes. As one possible simple test of the correctness of this 
algorithm, one can use it to reproduce the singular terms of the one loop helicity amplitudes 
of all 2 — > 2 subprocesses [0]. 

3 Singular terms of the one loop color subamplitudes 

In this section we give the singular soft and collinear terms of the one loop color subampli- 
tudes of the 2^3 subprocesses qQ — > qQg and qq — > ggg. One can obtain the result for the 
qq -^ qqg equal flavor subprocess from the unequal flavor process with simple permutations. 
We do not give the result for the subprocess gg -^ ggg since it was obtained earlier in 
and we found agreement. We note that we calculated the one loop helicity amplitudes of 
the subprocess qQ -^ qQg with standard diagrammatic method |^ and could confirm the 



correctness of the Born level algorithm for the singular terms also in this case. 

3.1 Subprocess class -^ qQQqg 

Here we present the singular terms of the one loop virtual corrections of the helicity color 
subamplitudes of the process —>■ qQQqg- The momenta of the partons are labeled as 

— ^ antiquarki (g) + antiquark2((5) + quark2(Q) + quarki(g) + gluon((7) . (3.1) 



The color structure of the amphtudes is the same at one loop as at tree level: 



A'''%q,Q;Q,q,g) = r 



7 . \^ )gi<|2'^g2giSig2(%; '^QJ '^Q; '^g; '^5) 

.(gi7^g2)e{'2,Q} 



(3.2) 
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The tree- level a\j subamplitudes are: 



0'if{hq,hQ;hQ,hq,+) = pa{hg,hQ;hQ,hq,+) 



m 



.(0) 



ttij {hg,hQ;hQ,hq,-) — pa{hq,hQ;hQ,hq,-) 






[w] [gj] ' 

The helicity dependence is completely absorbed in the factor p which is given by 

(mn)^ 



Pa{hq,hQ;hQ,hq,+) = Ti 
Pa{hq,hQ;hQ,hg,-) = ±i-. 



{qq){QQy 



(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 



'[qq][QQy 

where m and n labels those (anti) quarks which have opposite helicity to the gluon and the 
upper signs apply when m and n are both quarks or antiquarks, while the lower signs apply 
when one of them is a quark and the other is an antiquark. 

Using these expressions, the singular parts of the renormalized one-loop subamplitudes 
can be written in a simple form: 
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(3.9) 
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(3.10) 



The remaining subamplitudes can be obtained by simple label permutations 



a 



(l),sing (l),sin| 



qQ 



Qq 



a. 



(l),sing _ (l),sin: 



qq 



a; 



Ig^QjQ^g 



3.2 Subprocess class -^ qgggq 



(3.11) 



In this section we present the singular terms of the one-loop virtual corrections of the helicity 
color subamplitudes of the process -^ qgggq 
We use the the momentum labels as follows 

-^ quark(g) + gluon(l) + gluon(2) + gluon(3) + antiquark(g) . (3-12) 

The color structure of the amplitudes at tree level is given by 

(3.13) 



^^■"^^(g, g, g. g. g) = g' E (t^^t^^t^^),,^!"^^^, h, h, h„ h,), 

perm (1,2,3) 



while at one loop one finds 



A''''^{q,g,g,g,q) = / (£) 



E {T^'Ta'TS%,b'S3{hq,h,h2,h,,hg) (3.14) 

perm (1,2,3) 

+ E iT^%qSg,g,b^\K,h,h2,h„h,) 
i=l,2,3 

+ 5,,-Tr(T^^T^^T^«)6«(/i„ h, h, h, hg) 

+ 6ggTT{T^^TS'T^^)b'i\hg, h, h2, h, hg) 



In this equation k and / are the two indices from the set {1, 2, 3} \ i for a given i G {1, 2, 3}. 
The tree-level 6123 subamplitudes vanish for helicity configurations (±, hi, /i2, /is, ±) and 
(=F, h, h, h, ±). When there is one gluon with negative helicity (denoted by /) and two with 
positive hehcity, then 



r(0) 

^123l~5 'n, ''2, 



hi,h2,h3,+) 



bi23{+,hi,h2,h 



5 "-3) 



Pb{-,hi,h2,h3,+) 
(gl)(12)(23)(3g)(gg)' 

Pt{+,hi,h2,h,-) 
(gl)(12)(23)(3g)(gg)' 



(3.15) 
(3.16) 



where 



Pb{-,hi,h2,h3,+) 
Pbi+,hi,h2,h,-) 



-i{ql){ql)\ 



(3.17) 
(3.18) 



The subamplitudes for processes with opposite helicities can be obtained from ( |3.15| - |3.18|) 
by replacing ( ) with [ ]. 

Using these expressions, the singular parts of the renormalized one-loop subamplitudes 
can again be written in a simple form. For the ^L^ "^ amplitudes one finds 
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and the other permutations are obtained by simple rearrangement of indices: 



(3.19) 
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For the fo}^'''""'** amplitudes one finds 
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(3.24) 



(3.25) 



(3.26) 



4 Soft limits of the cross sections of 2 ^ 4 scattering 
processes 

The soft singular terms in the one loop radiative corrections of the 2 — > 3 hard processes 
determine the soft momentum dependence of a soft gluon emission from the hard partons 
that is the soft gluon momentum limit of general 2 — »• 4 processes 

parton(pi) + parton(p2) -^ parton(p3) + parton(j94) + parton(p5) + gluon(/c). (4.1) 



The soft momentum dependence factorizes into the eikonal factor 



e{m,n) = g 



2 Pn ■ Pr, 



Pn ' "^ Pm ' "^ 



SO that we may write in the soft limit |14, 15 



Y^ \M{pi,p2,P3,P4,P5,k)\'^^g'^ J2 e{m,n)ilj^n{pi,hi 

color 



(4.2) 



(4.3) 



n<m 



Integrating the (i-dimensional phase space over the gluon momenta we obtain the same 
singular factor as given by eq. (|2.7|) . Since the soft singularities cancel we can obtain 
the spin dependent '0nm(Pi, hi) coefficient functions of the eikonal factors by calculating the 
singular contributions of the virtual corrections to the colored summed cross section of the 
2^3 hard scattering process. 

On the other hand, the calculation of the spin dependent ipumiPi, hi) functions can also 
be carried out directly from the soft Bremsstrahlung equation (see ref. fl^). We carried 
out the calculation using both methods, in order to test the correctness of our results. It 
is an interesting feature of the soft limit that the spin dependence of the ipmn functions is 
factorizable: it can be absorbed into a simple factor p (see eqs. ( p.7| , |3.8| , p.l7| , p.l8] ) ) which 



is independent of the color structure. The hehcity dependence of the i'mniPi, hi) functions 
can be given as 

i^mniPt, hi) = \p{pi, /li)PVmn(Pi), (4-4) 

where ip^^ is the spin independent part of i'mniPi, hi). 

The angular ordering approximation to the soft gluon distribution in the hard process is 
derived by the azimuthal averaging over the eikonal factor. Shower Monte Carlo programs 
use this piece of information, therefore, it is of interest to know the soft coefficient functions. 
In particular one can develop an initial condition for the Monte Carlo programs where for a 
given 2 — i> 3 hard parton process a cone for the first gluon emission from each external line 
is chosen according to the probabilities defined with the help of soft emission. Therefore to 
obtain analytic expressions for the ipnm functions is of phenomenological interest. We note 
that 2-jet and 3-jet initial conditions are worked out for e~^e~ shower Monte Carlo programs 
but the 3-jet initial conditions are not yet developed for hadron-hadron and electron-proton 
collisions. 



4.1 Subprocess class -^ qQQqg{g) 

The momenta of the partons are labeled as in (^ 
the soft coefficients we obtain the expressions 
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(4.5) 



(4.6) 



(4.7) 



(4.8) 



Out of these four ip^^ functions one can obtain all the ten functions by changing the labels 
correspondingly 



^; 



qQ 



ip'tQ with q^ q,Q ^ Q 



• ^gb = ^Iq with q^Q,q^Q 

• i^fg = V'lg with q^q,Q ^Q 

• ^QQ = ^tq with q^Q,q^Q 

• V'gg = V"! with q^Q,Q ^q 

• V'Og = ^1 with q^Q,q^Q 

The spin summed coefficient functions are 

i^mn = —— (4, + 4, + 4,- + 4,-) Cn (4-9) 

4.2 Subprocess class ^ qQQQgig) 

In the case of equal flavor, we should use the same amplitudes as in the case of unequal flavor 
but we should antisymmetrize it in the momenta Q <-> g and in the corresponding color and 
helicity labels. 

The ipmn{hi) functions then can be written in the form 

-ipmnihi) = -j;^{\Pa{q,Q;Q,q;9)\'^ A^n + \pa{Q,q;Q,q;9)\'^ B^^ 



-2Re{pa{q, Q; Q, q; g)Pa{Q. q; Q, q; gyCmn)) • (4.10) 

To simplify the final expressions we use the notation 

m - ^ f, = ^^ 



S^gSjg {W){9J) 

but here g denotes the hard gluon momentum, and i,j can run over values q,Q,q or Q. 



Then the auxiliary functions A, B, C appearing in ( [4. 101) can be written in the simple form 
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9 



^95 = ~j;p [f'QQfQq + fQQfqQj + Jp yfqqfhq + fmfqQ 

~ fqqfqq + ^ fqQfqQ ~ ^c [fqQfgQ + fQqfQQ) 

+ yfqQJQq + JQqJQQ + JqQJQQ + JQqfqq + JqQJqq + /q?/, 

-4,-0 = (^ - y) (W) + {«« - {«Q} 
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c c 

~ fQqfm + ^fqOfqQ ~ ^c [fqqfqq + fqQfQQj 

+ [fQqfQQ + fqqfqQ + fqQfqQ + /Qg/gg + fQQfqq + fqQfqq) 
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^95 ~ ^c yfQqfqq + fqqfqq + fqqfqq fqqfqq fqqfqq fqqfqq 
"•" Uqqjqq JqqJqq JqqJqq JqqJqq 



In order to get the expressions for the other ?/^^„ functions, we have to do the same relabehng 
procedure as in the previous process now with the A, B and C functions . The spin summed 
V'nin functions can be obtained by 

\Pa{q,Q-,Q,q-,g)? -^ (^'q + 4q + ^'q + 4q) (4-11) 

SqqSqq 
2 

|Pa(Q,g;Q,g;^)r ^ (^q + 4q + 4? + ^q) 

^qQ^qq 

Re(pa(g, Q; Q, ?; 9)PaiQ, q; Q, q; 5')*C'mn) 

g"^ \- s"^ 

~^ T~o I I {SqqSQq + SqQSgQ - SqQSqQ)Re{C^ri) 
^qq^qq^qq^qq 

4.3 Subprocess class -^ qqgggig) 

We use the the momentum labels as in (p.l2|) and we introduce the short hand notation 

"P = Si2Si3S23SiqSiqS2qS2qS3qS3qSqq. 



10 



Then for the soft color coefficient functions we obtain 

Ci = — — (4.12) 

Nc SiqSiqS2qS2qS3gSsg 



S23SlgSgg \S12S3g Si3S2qJ 

NcV r 2 2 ^2 

~ jy |~<J ■ S23SigS2qS3q — S23SiqSiqSqq + S23S^^S2qS3q 

+ [si3SlqS2qS2qS3q + S23SlqSiqS2qS3q + Si3S23SlqS2qSqq 

+ sl^S2qS2qSqg + Si3Sigsl^S3q - Si3SiqS3qsl^ +(2^3) | . 

From the above expression we get five other ip^^ functions by changing the labels of the 
momenta in the spin independent part 



. ijf^ ■.{q,q, 1,2,3)- 


^{q,q,2,3,l) 




• ^IJ%: iq,q, 1,2,3)- 


^iq,q,3,l,2) 




• ^|\ :(g,g, 1,2,3) - 


^{q,q, 1,2, 3) 




• ^1 : i<l, ^, 1, 2, 3) - 


^{q,q,2,3,l) 




. ^|'3:(g,g,l,2,3)- 


^{q,q,3,l,2) 




There are still two additional independent coefficient functions which cannot be 
symmetry properties 


obtained by 


fqq 


(1 - N^) 




N^ SlgSlgS2gS2qS3qS3q 




+ 


N^V „ r 2 2 

> , 2 ■ Si2S^gSlqS2q 2 ■ Sl2S^gSlgS2q 

^ Z(123) 
+ \S\2S3qS3q + S\3S23Sqq\ \S\qS2q + S\qS2q\ 


(4.13) 


— 


^ ^ / , S\2S\3S2qS3qSqq 2 ■ S23S\qS\qS2qS3q 
^ ^^c Z(123) 





- 2 ■ s\gS23S2qS3q + Sl2Sl3S23SlgSlg + (^ ^ ?) 

where Z(123) denotes the cyclic permutations and 

N'^V ( 1 1 \ 

^l2 = ^^ + + {q--q)] (4.14) 

Sgg \-Sl2Sl3-S2qS3q Sl2S23-^lq-S3q ) 

- — ^{-2 ■ S12S13S23 [SlgSig + S2<?S2gJ 

—2 ■ SY2^3qS3qSgg — Si3S23Sgg \S\gS2q + S\qS2q\ 
+2 ■ S23 {S\q + S15] {SxqS2qS3q + SlgS2gS3g] 
+ 2 ■ Si3 {S2q + S2g] {s\qS2qS3q + SlgS2gS3g]} 

^mn = 2 5] (4S^5 + S.,4)C, (4.15) 

JG{1,2,3} 
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Spin summation results in 



4.4 Subprocess class -^ 99999(9) 

We introduce the simple notation for labeling the quarks and gluons as 

-^ gluon(l) + gluon(2) + gluon(3) + gluon(4) + gluon(5) (4.16) 

Due to Bose symmetry it is sufficient to calculate only one of the ten ipmn functions 

i^fs = 2iV^V Y: ^ (4-17) 

P(2,4,5) '5l3"S32"S24'S45'S51 



.S12S14S15S24S25 S14S15S23S25S34 
1 1 

+ 



+ 



S14S15S23S24S35 S12S14S15S34S35 

1 1 

+ 

S12S24S25S34S35 S12S15S23S34S45 

1 2S13 



S12S14S23S35S45 S12S14S15S23S34S35 
S24 S25 



■Sl2'5l4'523'525'S34S45 ■S12S15S23S24S35S45 

Although in the expression for ^^3 the symmetry 1 ^^ 3 is not manifest, using momentum 
conservation we can check that it has the symmetry 1 ^^ 3. We get the other ?/'^^ functions 
now simply by the usual replacements of the labels of the momenta, i.e. ipf^ -^ ijjf2 by 3 ^^ 2. 
The spin dependent part of the born amplitude is given by 

Pd = i(mn)^ 

if m and n are the only negative helicity gluons and by 

Pd = i[mn]^ 

if m and n are the only positive helicity gluons. So, for the helicity summed ipmn functions 
we get 

^P^n = 2Y.st^rL■ (4.18) 

i<j 

5 Conclusions 

We presented an algorithm of general validity, which enables us to calculate the singular 
terms of any one-loop helicity amplitude in QCD. The correctness of the algorithm was 
tested by explicitly demonstrating the cancellation of the soft singularities between the loop 
and Bremsstrahlung contributions for spin dependent 2 — >■ 3, 2 -h> 4 processes. We also 



emphasized the universality of the collinear singularities (see (2^)). Furtheron, we obtained 
short formulae for the soft limit of the spin dependent cross sections of all 2 — * 4 processes. 
These results represent a necessary input for shower Monte Carlo programs for the production 
of three well-separated jets. 
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